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(57) Abstract 



Method for the measurement of the cardiac flow, in which the volume of blood ejected from the left ventricle (LSV) or the volume 
of blood ejected from the right ventricle (RSV) is expressed and calculated as a function of at least one of the independent contributing 
factors of the area under the pressure curve recorded by means of an appropriate sensor and as a function of the hydraulic impedance, and 
in which the cardiac flow Q is calculated according to the relation Q - LSV (or RSV) * HR, where HR is the heart rate* 
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METHOD AND APPARATUS FOR MEASURING CARDIAC FLOW OUTPUT 
Field of the invention 

The present invention refers to a method and an apparatus for determining the 
stroke volume - i.e., the volume of blood expelled from the left ventricle (LSV), the 
volume of blood expelled from the right ventricle (RSV) - and hence the cardiac 
output Q - i.e., the stroke volume multiplied by the heart rate (HR) in a 
continuous way, using invasive and non-invasive indirect techniques, so as to 
enable acquisition of this important haemodynamic parameter in various clinical 
and non-clinical situations, as well as in the course of ergometric tests. 
Prior art 

For the measurement of cardiac flow Q, the invasive methods that are currently 
most extensively used are the Thermodilution Method (TDM), Fick's Method (FM), 
and a method that uses the arterial pressure signal p(t) measured in the aorta or in 
the pulmonary artery, referred to as the Pulse Contour Method (PCM). 
This method which uses the signal pressure is not very reliable and for this reason 
it is necessary to make a calibration. This is usually the TDM. At the present time 
reliable results have not been attained by this method. 

This PCM method derives from an original idea of Herd [Herd J.A. et a/., 1864] and 
from the theory referred to as the Windkassel (German for "air chamber") theory of 
Franck (Franck O., 1930), and is based on the existence of a relationship between 
the volume of blood expelled by the left ventricle (LSV) or the volume of blood 
expelled by the right ventricle (RSV ) f and the area under the pressure curve p(t). 
The fundamental relation used for calculating the stroke volume is SV=A/Z0 where 
A, expressed in [mmHg t], is the area under the pressure curve p(t) (see figure 
A1), and Z0, expressed in [mmHg/cm/t], is the hydraulic impedance which 
depends upon the dynamic resistances and upon the compliance of the artery 
wall. LSV is measured in [cm 3 ] (see figure A1), hence Q = LSV*HR is the cardiac 
flow expressed in litres per minute if the heart rate is measured in beats per 
minute. In this connection we recall that the plot of the arterial pressure with 
respect to time is determined by the magnitude of LSV and by the vascular 
impedance. Consequently, the Pulse Contour Method endeavours to separate and 
analyse these two contributions; however, the method is unable to determine the 



WO 00/64339 



PCT/EPOO/03697 



2 

two contributions as independent variables overtime. 

Applying the Windkassel theory, many studies have attempted to determine LSV 
only from the pressure waveform and from the characteristics linked to 
transmission of the wave in the aorta or in the pulmonary artery [Remington J.W. 
et a/., 1948; Warner H.R. et a/., 1953; Herd J.A. ef a/., 1966; Kouchoukos N.T. et 
a/., 1970], 

The original idea of Franck has subsequently been applied over the years and has 
made it possible to estimate LSV in a continuous way from the measurement of 
the pressure signal in the aorta or in the pulmonary artery [McDonald D.A. et a/., 
1974; Wesseling K.H. ef a/., 1976; Tajimi T. et a/., 1983; Wesseling K.H. et a/., 
1993]. 

However, in the concrete application to the various possible clinical situations, the 
Pulse Contour Method currently requires a "calibration" for calculation of the 
hydraulic impedance. For the calibration, generally one of the other two methods 
referred to above, i.e., the thermodilution method and Fick's method, is used, or 
else linear regressions of aortic parameters, such as the diameter of the aorta, and 
the age, sex, height and weight of the patient are used. 

Unfortunately, the calibration and regression factors are subject to error, given that 
the methods from which they are obtained are in turn imprecise and that the 
regressions are in any case obtained on a limited number of subjects, and are 
hence acceptable only as a mean and not as a true measurement of the quantity 
investigated. 

In fact, cardiac flow estimated using the thermodilution method and Fick's method 
are not always in agreement with the clinical parameters obtained using other 
diagnostic techniques, and this mainly occurs in patients suffering from certain 
forms of heart disease, such as disease involving dilatation of the heart, valvular 
cardiopathy, and cardiac fibrillation. 

To provide an example, consider two possible signals in the aorta studied between 
. the points of opening and closing of the ventricle. These signals will in general 
present the same area but different forms, with different times of attainment of the 
systolic point. 

The traditional Pulse Contour Method will therefore yield the same exact 
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measurement (same integral) evaluated on the basis of the calibration impedance. 
However, it is evident that from a different form of the signal there must derive a 
different impedance, which cannot be evaluated. 

Hence, the limits of the invasive techniques currently in use are: a) the poor 
5 precision achievable in the estimation of cardiac flow on account of clinical 

illnesses; b) the non-applicability in general on account of the pathological 

condition of the patient; and c) the impossibility of applying the said invasive 

techniques, for example during an ergometric test. 

Scope of the invention 
10 A first scope of the invention is to obtain measurements in a continuous way that 

are more reliable than the ones obtained using the invasive and non-invasive 

techniques currently applied. 

A second scope is to render the measurement substantially independent of the 
point of application of the sensor by introducing variations in the specific formulas, 
is without any need for prior calibration of the measurement. 
Summary of the invention 

The above purposes have been achieved according to the invention using a 
method which directly obtain the cardiac flow from the pressure signal measured 
in an invasive way, in the ascending aorta, in the pulmonary artery and in femoral 

20 brachial and radial, or measured in a non-invasive way, for example from the 
arteriole of the finger using a cuff meter. According to the method, the impedance 
of the pressure signal is calculated on the basis of the points of resonance of the 
signal by assimilating the signal to that of a flow in an elastic tube and assuming 
constancy of Young's modulus, which is in particular taken to have a value of 

25 unity. In this way, it is possible to calculate the cardiac flow without any longer 
having to use various calibrations, but exclusively from an analysis of the pressure 
wave and its characteristics. 

Preferably, the hydraulic impedance is calculated by means of an analysis of the 
first and second derivatives with respect to time of the pressure signal recorded. 
30 According to a further aspect of the invention, a correction is moreover made of 
the mean pressure value to be used for the calculation of LSV in order to take into 
account the attenuations of the said value in the various points of possible 
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recording of the signal. 

According to a further aspect of the invention, from the signal recorded in the 
finger (or from some other point in a non-invasive way), the method makes 
possible a direct reconstruction of the signal in the aorta and in the pulmonary 
artery, and from the latter signal a reconstruction of the cardiac flow. 
More specifically, according to the invention In order to obtain the SV estimation 
we have taken into consideration the wave pressure in the ascending aorta and/or 
in pulmonary artery, the artery compliance (E) and the periferical resistance (R). 
Therefore we have taken into consideration that 1) the SV depends on the 
pressure variation obtained at the opening of the ventricular valve (that is the 
difference between systolic and diastolic pressure divided by the time necessary to 
pass between the systolic and diastolic) and 2) the SV is conditioned by E and R. 
In order to obtain these contributions we have taken into account the value of the 
dicrotic pressure and the other characteristic points between the systolic dicrotic 
pressure (this pressure values must be divided by time. This time is the difference 
between the end of the cardiac beat and the time of the event being taken into 
consideration). 

Therefore we have considered SV as being determined by three points: 1) the 
bolus of blood ejected by the ventricular; 2) the reaction due to the aortic wall; 3) 
resistance due to periferic arterial cycle. As the value of the pressure at the point 
where it is taken is the result of these three components at the same time we have 
studied our system in a perturbative way. Therefore we have considered the 
principle contribution of the ventricular and of the system E and R, the first being 
given by 1) as described above and the second, the E and R system principally 
contributes to the closure of the valve ( point of diacrotic pressure) This last event 
point is conditioned by a series of perturbations on the pressure signal after the 
cardiac valve, according to the vase being crossed and the length of the course. 
That is it is necessary to take into consideration not only the contribution due to 
the systolic and to the diacrotic above described, and when present those due to 
secondary perturbations. 

In conclusion all the event points which have been taken into consideration are the 
moments in which there is a state of balance between the various points (blood 
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ejected from the ventricular-E-R): the "principle" balancing points ( systolic and 
dicrotic points) can be or not "accompanied" by other balancing points (how and if 
to analyse them will be described later). All this information can be found in the 
wave pressure which flows after having been generated by the ventricular (right or 
left). 

Advantageously, with the method according to the invention it is possible to 
establish a relationship between the hydraulic impedance and the usable time, 
also in combination with known methods (e.g., the thermodilution method) which 
involve a phase of calibration of the signal recorded, where the contribution of the 
area under the pressure curve is considered variable over time, whilst the 
contribution of the impedance can only be considered as constant. 
In particular, by the method of the present invention (hereinafter called pulse 
analytical method (PAM)) it has been possible: a) to find the SV from the signal 
pressure recorded invasively in ascending aorta and in pulmonary artery; b) to find 
the SV from the arterial signal pressure invasively recorded (brachial, radial, and 
femoral artery) and non invasively recorded (e.g. from the pressure obtained with 
the oscillometry method from the arterial finger). 

In this way we estimated the LSV and RSV and so we determined the true Q in a 
way that is completely free from any calibration. Therefore these results are 
obtained only by the analysis of the wave pressure ( depending only on where the 
wave pressure was taken). 

According to the invention an apparatus able to perform the method is provided. 

The apparatus comprises a microprocessor unit able to receive the blood pressure 

signal and to analyse it over the time in order to determine the above identified 

parameters and to calculate the cardiac output Q therefrom. 

In a preferred embodiment the apparatus further comprises a sensor in the shape 

of a cuff meter able to be applied to a finger and to acquire the blood pressure 

signal. 

Brief description of the figures 

- Figure A1 shows a cardiac pressure signal diagram as analysed by the prior art 
method; 

- Figures 1-19 show the form of the cardiac pressure signal and of its first and 
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second derivatives, the signal being taken in various recording points; 
- Figure 20 shows a reconstruction, according to the present method, of the signal 
in the aorta starting from the pressure signal recorded at the arteriole of the finger. 
Detailed description of the invention 
5 With reference to the annexed figures, various examples of application of the 
method are described in what follows. 
Example I 

Aj Relation between LSV and pressure jaken in the ascending aorta (Pulse 
Analytir^l Method. Ao rtic: PAMA1 (Fia.1-6^ 
io i) The PAMA determines the heart flow Q in litres per minute using the following 
general relation (the pressure signal is acquired in the ascending aorta at 1000 

Hz): 

LSV = [K[A/((Za1+Za2)*1000)+A/((Za1+Za2)*1000r(Pm-K1)/K1]]/1000 
Eq. [1] 

15 where: 

K = 1 and has the dimensions [(\m*sqrt(2p/(p)*Vm], expressed in [I 3 / 1 2 ]; 
Xm is the mean wavelength, approximately equal to 10 m 
Vm is the mean velocity, approximately equal to 10 m/s 
p is the blood density; 

20 A is the integral between t1 (time at the diastolic pulse in [ms]) and tdic (time at the 
dicrotic pressure in [ms]) under the pressure curve p(t), expressed in [mmHg*ms] 
(Figure 1); 

K1 = 100, expressed in [mmHg], and represents the correction factor of the mean 
pressure; 

25 Za1 = (psys-p(1))/tsys, expressed in [mmHg/ms]; 
Za2 = (pdic/tfinal - tdic), expressed in [mmHg/ms]; and 
Pm = (psys+2p(1))/3. see the following Notel 

tfinaN time of the beat being taken into consideration (time origin in t1 and final in 
tfmal) 

30 The cardiac flow is thus 
Q = LSV*HR 

where Q is expressed in [lit/min]; 
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HR = 60000/T; and 

T is the cardiac period expressed in [ms]. 

This relation was applied in cases in which the pressure curve and the 
corresponding mean of the tangent at 21 points (i.e., the first derivative d") and the 
mean of the tangent at 21 points of the mean tangent (i.e., the second derivative 
d") were those shown in Figures 2 and 3 and could be associated to the recording 
point. 

in With -Za3 

In the cases where the pressure curves in the ascending aorta were of the type 
shown in Figure 4, and the corresponding first and second derivatives, d' and d". 
were as those shown in Figures 5 and 6 and revealed the point of resonance at 
time t3, then the relation became: 

LSV = [K[A/((Za1+Za2-Za3)*1000)+A/((Za1+Za2-Za3)*1000r(Pm-K1)/K1]] /1000 
Eq.[2] 

where the symbols have the same meaning as in Eq. [1], and where t3 is the time, 

in [ms], at the minimum value of d" between the time tsys and the time tdic, and p3 

is the corresponding pressure expressed in [mmHg] at time t3 (see Figure 6) and 

Zf3 = (P3/(tfinal -t3) mmHg/ms 

In a similar way it is possible to calculate Q = LSV*HR. 

Note 1 

The mean pressure for the pressure measured in the ascending aorta must be 
considered as such for the interval 90-110 mmHg; for the mean pressure 
between 1 10 - 120 and 90 - 80 mmHg it must be considered at 50% (for example 
for a Pm= 118 mmHg for our method is = 114 mmHg); for mean pressure values 
between 120- 130 and 80 - 70 mmHg it must be considered at 25%, for mean 
pressure values>=130 and <=70 mmHg it must be considered 13% 
Example II 

Relation between RSV and the pre ssure taken in the pulmonary artery (Pulse 
Analytical Method. Pulm onary: PAMP) 

Relationship between the volume of blood expelled from the right ventricle RSV 
and the pressure measured in the pulmonary artery. The corresponding signal 
pressure is similar to the one represented for aortic pressure, but for variations in 
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scale (see Figure 7). 

The PAMP determines the heart flow Q in litres per minute using the following 
general relation (the pressure is acquired in the pulmonary artery at 1000 Hz): 
i) Case with mean pressure in pulmonary artery >= 19mmHg 
5 RSV = [K[A/((Za1+Za2ri000)+A/((Za1 + Za2ri000)*(Pm-K1)/K1]]/1000 [Eq.3] 

where: 

K = 1 and has the dimensions [(X m *sqrt(2p/(p)*V m ], expressed in [I 3 / 1 2 ], p being the 
density of the blood; — 

A is the integral between t1 (time at the diastolic pulse in [ms]) and tdic (time at the 
10 second dilatation of the artery in a dicrotic pulse in [ms]) under the pressure curve 

p(t), expressed in [mmHg*ms]; 

K1 = 12, expressed in [mmHg]; 

Za1 = (psys)/tsys, expressed in [mmHg/ms]; 

Za2 = (pdic/tfinal-tdic), expressed in [mmHg/ms]; and 
is Pm = (psys+2p(1 )/3); see following Note 2 

Q = RSVHR, 

where Q is expressed in [lit/min]; 
HR = 60000/T; and 

T is the cardiac period expressed in [ms]. 
20 In Figure 7 a signal acquisition of the pressure in the pulmonary artery is shown. 
For the pressure in pulmonary artery we have variations for d' and d" similar to 
those of the aorta. In consequence the determination of the point of dicrotic 
pressure (Pdic), the systolic pressure(Psys), diastolic pressure (P(1)) and the 
relative times is as described above. 

25 ii) Case with mean pressure in pulmonary artery <= 19mmHg 
In the cases in which Pm is <= 19 mmHg the relation becomes: 
RSV=[K [ A / ((Za1+Za2)*1000)+A / ((Za1+Za2)*1000]] / 1000 Eq.[4] 
With the same meaning for the symbols as in the preceding cases. In the same 
way it is possible to calculate Q=RSV*HR 

30 Note 2 

The mean pressure in the case of pressure taken in the pulmonart artery must be 
taken as such for the interval of pressure between 19-28 mmHg; for values of 
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mean pressure between 28 - 33 mmHg it must be considered at 50%, for values 
of mean pressure > 33mmHg it must be considered at 25% (for example a 
pm=43mmHg for our method is equal to 33 mmHg); for values < 19mmHg we are 
in case ii) and therefore we do not use the mean pressure. 
Example III 

CV Relation between LSV and the pressure non invasiv elv recorded from the 
arterial finger (Pulse Analysis Metho d Finger: PAMR 

Direct relationship 

i) The PAMF determines the cardiac flow Q in litres per minute using the following 
general relation (the pressure is acquired at the finger of the left hand at 1000 Hz): 
LSV = [k[A/((Zf1+Zf2)*1000)+A/((Zf1+Zf2ri000r(Pm-K1)/K1]]/1000 Eq. [5] 

where (Figure 8): 

K = 1 and has the dimensions [(Xm*sqrt(2p/(p)*Vm], expressed in [I 3 / 1 2 ]; 
A is the integral between t1 (time at the diastolic pulse in [ms]) and tdic (time at the 
dicrotic pressure in [ms]) under the pressure curve p(t). expressed 
in [mmHg*ms]; 

K1 = 90, expressed in [mmHg]; 

Zf1 = (psys-p(1))/tsys, expressed in [mmHg/ms]; 

Zf2 = pdic/(tfinal - tdic), expressed in [mmHg/ms]; and 

Pm = (psys+2p( 1 ))/3 . see following Note 3 

The corrected volume of blood expelled from the left ventricle (LSVC) is 

LSVC = [LSV+LSV*abs(delta(Pd1-pdic))/(psys-pdias)] [6] 

where: 

(Pd1-pdic) = variation of pressure of the dicrotic point (Pdic) at its maximum 
(Pd1)=[mmHg]. This correction exists only when there is an increase in the 
pressure after the dicrotic pressure ((Pd1 - Pdic)>0). In the cases in which the 
increase in pressure is not present ((Pd1 - Pdic)<=0) we have LSV=LSVC. 
Psys is the systolic pressure, expressed in [mmHg]; 
Pdias is the diastolic pressure, expressed in [mmHg]; and 
the term Pd1 is calculated immediately after the dicrotic point and is the maximum 
value of the curve after (Pdic). 
Q = LSVC*HR 
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where Q is expressed in [lit/min]; 
HR = 60000/T; and 

T is the cardiac period expressed in [ms]. 

The above relation was applied in the cases where the pressure curve and the 
5 corresponding first and second derivatives, d' and d\ were those shown in Figures 
9 and 10. 
ii) With -Zf3 

In the cases where the pressure curves were of the type shown in Figure 1 1 and 
the conesponding first and second derivatives, d' and d", were as those shown in 
10 Figures 12 and 13, the relation became: 

LSV = [k[A/((Zf 1 +Zf2-Zf3)*1 000)+A/((Zf 1 +Zf2-Zf3)* 1 000)*(Pm-K1 )/K1]]/1 000 

Eq.[7] 

where: 

Zf3 = P3/(tfmal - 13); and 
15 the symbols have the same meanings as specified previously, and t3 is the time, 

in [ms], of the minimum value of d" between the time tsys and the time tdic, and P3 

is the corresponding pressure, expressed in [mmHg] at time t3 (see Figure 11). 

LSVC = LSV+LSV*abs ((Pd1-Pdic))/(psys-P(1 )) Eq. [8] 

In a similar way it is possible to calculate Q = LSVC*HR. 
20 iii) With -2Zf3 

In the cases where the pressure curves were of the type shown in Figure 14 and 

the corresponding first and second derivatives, d' and d", were as those shown in 

Figures 15 and 16, the relation became: 

LSV = [k[A/((Zf1 +Zf2-2Zf3)*1000)+A/((Zf1 +Zf2-2Zf3)*1 000)*(Pm-K1 )/K1]]/1 000 

25 Eq.[9] 

where Zf3 = P3/(tfinal - 13); and 

the symbols have the same meanings as specified previously, and t3 is the time, 
expressed in [ms], of the minimum of d" between the time tsys and the time tdic, 
and P3 is the corresponding pressure, expressed in [mmHg] at time t3 (see Figure 

30 14). 

LSVC = LSV+LSV*abs ((Pd1-Pdic))/(Psys-P1) [10] 

In a similar way it is possible to calculate Q = LSVC*HR, expressed in litres per 
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minute. 

iv) With -2Zf3-Zf5 

In the cases where the pressure curves were of the type shown in Figure 17 and 
the corresponding first and second derivatives, d' and d", were as those shown in 
5 Figures 18 and 19, the relation became: 

LSV = [k[A/((Zf1+Zf2-2Zf3-Zf5)*1000)+A/((Zf1 + Zf2-2Zf3- 
Zf5)*1 000)*(Pm-K1 )/K1]]/1000 

where Zf3 = P3/(tfinal - 13) — - 
where Zf5 = P5/(tfinale-t5) 
10 the symbols have the same meanings as specified previously, and t5 is the time, 
expressed in [ms], of the minimum of d" between the time tsys and the time tdic, 
and P5 is the corresponding pressure, expressed in [mmHg] at time t5 (see Figure 
17). 

LSVC = LSV+LSV*abs ((P1-Pdic))/(Psys-P1)[11] ex 10 
15 In a similar way it is possible to calculate Q = LSVC*HR, expressed in litres per 
minute. 
Note 3 

The mean pressure in the case of the pressure taken at the arterial finger non 
invasively must be considered as such for the interval of mean pressure between 

20 70 - 1 10, for the values of mean pressure between 1 10 - 1 50 and 70 - 40 mmHg 
it must be considered at 50% (for example a pm=128 for our method is = 119 
mmHg); for mean values of pressure >150 and < 40 mmHg it must be considered at 25% 

Reconstruction of the pressure signal in the ascend ing aorta bv means of linear 
multiple regression in the time domain, bv use of Zf1-Zf5 

25 For these reconstruction, basically linear multiple regressions were used. In order 
to reconstruct the signal recorded in the ascending aorta (or in the pulmonary 
artery) using a cardiac catheter from the arterial signal recorded in a continuous 
way by means of a small cuff wrapped around the middle finger of the left hand, a 
linear multiple regression was used in which the reconstructed pressure signal 

30 was obtained in two successive steps: 

1) An estimate was made of the mean pressure during the cardiac cycle in the 
ascending aorta (or in the pulmonary artery) from the signal taken at the finger, 
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deriving the value Pmf (mean pressure in the aorta estimated from the recording 
taken at the finger) from the formulas used in the various cases of analysis of the 
arterial signal referred to in the previous points: 
Pmf = LSV*Ztot/(k*A) [11] 

2) The waveform in the ascending aorta (or in the pulmonary artery) was 
reconstructed from a fit that used the following parameters: 
y = aO*Pmf + a1*fin + a2*abs(derfin) + a3*abs(der2fin) + a4*abs(der3fin) + 
a5*(intfin) + a6*slope*abs(derfin) + a7*slope*zZf1 + a8*slope 
+ a9*maxfin + a10*minfin + a1 1*HR*(intfm(up to the point considered))* 
a12*areaf + a13*zZf1 + a14*zZf2 + a15*zz3f + a16*zz4f a17*Zf5 
areaf = cof*(Zf1 +Zf2) / (Zf1 +Zf2-n*z3f- Zf5) [12] 
where 

Zf5 and n = 0, 1 and 2 according to the criteria described previously; 
zz4f = Pd1/(tfinale-td1 ) (Figure 14); 

- fin is the pressure at the finger; 

- abs(derfin) is the absolute value of the first derivative in the pressure point 
considered; 

- abs(der2fm) is the absolute value of the second derivative in the pressure point 
considered; 

- abs(der3fin) is the absolute value of the third derivative in the pressure point 
considered; 

- infin is the integral up to the point considered of the signal at the finger; 

- slope is the angle between the horizontal axis and the straight line passing 
through the minimum points on the left and on the right of the cardiac cycle; 

25 - maxfin and minfin coincide with the systolic pressure and the diastolic pressure; 

- areaf is the total area of the pressure signal; 

and the remaining symbols have the same meanings as specified previously. 
A number of reconstructions obtained are illustrated in Figure 20. 
The errors in the comparison between the reconstructed curve of the signal 
30 registered non invasively and that taken directly near the ascending aortic are 

SD(mmHg) Max(mmHg) Min(mmHg) 

1.16+5.67 2.38-16.40 -2.82 + -16.41 
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mean: 3.41 9.37 -9.32 

With SD= Standard Deviation: the minimum of the interval is obtained for the 
riconstrruction of the points around the diastolic pressure, the maximum of the 
difference is obtained near the point of the systolic pressure. 
Max= interval of over estimation of the pressure in the point taken into 
consideration reconstructed and that actually measured with the catheter during 
the cardiac beat: the minimum of this interval is obtained for the reconstruction of 
the points around the diastolic pressure, the maximum of the difference is obtained 
near the points of systolic pressure. 

Max= interval of over estimation of the pressure in the point taken into 
consideration reconstructed and that actually measured with the catheter during 
the cardiac beat: the minimum of this interval is obtained for the reconstruction of 
the points around the diastolic pressure, the maximum of the difference is obtained 
near the points of systolic pressure. 

Min= interval of under estimation of the pressure in the point taken into 
consideration reconstructed and that actually measured with the catheter during 
the cardiac beat: the minimum of this interval is obtained for the reconstruction of 
the points around the diastolic pressure, the maximum of the difference is obtained 
near the points of systolic pressure. 

What is important in this calculation are Zf 1 , Zf2,Zf3, Zf5. which we considered in 
point C): these are necessary to have satisfactory results. 
D) Relation between LSV and the pressure recorded invasively from femoral artery 
or from another periferical point such as brachial or radial artery (Pulse Analytical 
Method , Brachial, Radial and Femoral, PAM(BRF) ).... 

For these case we have seen that the formula to use are of the type used in the 
case C) non invasively with the following precisions: i) K1 for these invasive 
signals must be considered =100; ii) note 3 remains unchanged. 
According to the invention, the method can be applied in combination with known 
methods (such as the thermodilution method) comprising a phase of calibration of 
the signal recorded, in which the contribution of the area under the pressure curve 
is considered variable over time, whereas the contribution of the impedance can 
only be considered constant. 
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In this case, the proposed method also makes it possible to take into account even 
major variations in the heart rate, in the pressure values and in the pressure 
waveform for purposes of calculation of the impedance. 

It may therefore be concluded that, both in the case of normal subjects and in the 
5 case of patients affected by various pathological conditions, the method proposed 
represents an effective and advantageous diagnostic tool in the invasive and non- 
invasive evaluation of cardiac flow. 

In addition, the method can be applied both in healthy subjects and in subjects 
presenting cardiocirculatory alterations who are undergoing ergometric tests that 

10 are aimed at establishing the level of haemodynamic response to the tests. 

It should be emphasized that the present method is based only on the study of the 
pressure signal (recorded invasively in the pulmonary artery and in the aortic arch, 
or in any other major arterial vessel, or else non-invasively at the finger), and is 
independent of the anthropometric data and of the age of the subject examined. 

15 The present invention further comprises an apparatus for measuring cardiac flow, 
comprising at least one sensor for detecting a blood pressure signal and a 
computer unit connected to the said sensor for the execution of a measurement 
according to the above described method and provided with at least one output 
device of the measured value. 

20 Preferably, the apparatus comprises a storage medium containing a computer 
program to execute a method according to at least one of claims 1 to12. 
The invention further relates to a computer program loadable in a computer unit in 
order to execute the method. 
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CLAIMS 

r 1 . Method for the measurement of the cardiac flow output, comprising the steps of: 

2 - recording of the blood pressure signal by means of an appropriate sensor; 

3 - calculating over the time the volume of blood ejected from left ventricle (LSV) or 

4 the volume of blood ejected from the right ventricle (RSV) as a function of the 

5 artery compliance and/or of the periferical resistance and/or of the hydraulic 

6 impedance; 

7 - calculating the cardiac flow output (Q) according to the relation Q = LSV (or RSV) 

8 * HR, where HR is the heart rate. 

1 2. Method according to Claim 1, in which for the calculation of LSV a correction of 

2 the mean value of the pressure is performed according to the site chosen for the 

3 recording obtaining a corrected value LSVC. 

1 3. Method according to Claim 1 , in which the analysis of the blood pressure signal 

2 envisages the study of the first and second derivatives of the signal for calculating 

3 the resonance points from which the hydraulic impedance may be calculated. 

1 4. Method according to Claim 1 , in which the pressure is acquired in the ascending 

2 aorta, and LSV is calculated by (Eq. [1]). 

1 5. Method according to Claim 1 , in which the pressure is acquired in the ascending 

2 aorta, and LSV is calculated by (Eq. [2]) 

1 6. Method according to Claim 1, in which the pressure is acquired in the 

2 pulmonary artery, and RSV is calculated by (Eq. [3]). 

1 7. Method according to Claim 1, in which the pressure is acquired in the 

2 pulmonary artery, and RSV is calculated by (Eq. [4]). 

1 8. Method according to Claim 2, in which the pressure is acquired in a non- 

2 invasive way from the arteriole of the finger, and LSVC is calculated by (Eq. [5, 

3 6]). 

1 9. Method according to Claim 2, in which the pressure is acquired in a non- 

2 invasive way from the arteriole of the finger, and LSVC is calculated by (Eq. 

3 [7,8]). 

1 10. Method according to Claim 2, in which the pressure is acquired in a non- 

2 invasive way from the arteriole of the finger, and LSVC is calculated by (Eq. 

3 [9,10]). 
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1 11. Method according to at least one of the foregoing claims, in which the pressure 

2 is acquired in a non-invasive way from the arteriole of the finger, and a 

3 reconstruction is carried out of the pressure signal in the ascending aorta by 

4 means of a linear multiple regression. 

1 12. Method according to Claim 11, in which the said reconstruction consists in an 

2 estimation of the mean pressure during the cardiac cycle in the ascending aorta 

3 from the signal recorded at the finger, deriving therefrom the value 

4 Pmf = LSV*Ztot/(k*A)(Eq.[11]) 

5 and in the reconstruction of the waveform in the ascending aorta by means of a fit 

6 which uses the following parameters: 

7 y = aO*Pmf + a1 *fm + A2*abs(derfin) + a3*abs(der2fin) + a4*abs(der3fin) + 

8 a5*(intfin) + a6*slope*abs(derfin) + a7*slope*zZfl + a8*slope 

9 + a9*maxfin + a10*minfin + a1 1*HR*(intfin(up to the point considered))* 

10 a12*areaf + a13*zZf1 + a14*zZf2 + a15*zz3f + a16*zz4fa17*Zf5 
n areaf = cof*(Zf1+Zf2)/(Zf1+Zf2-n*Zf3-Zf5) (Eq. [12]) 

12 where 

13 Zf5 and n = 0, 1 and 2; 

14 zz4f = Pd1/(tfinal-td1); 

15 fin is the pressure at the finger; 

16 abs(derfin) is the absolute value of the first derivative in the pressure point 
n considered; 

18 abs(der2fin) is the absolute value of the second derivative in the pressure point 

19 considered; 

20 abs(der3fin) is the absolute value of the third derivative in the pressure point 

21 considered; 

22 infin is the integral up to the point considered of the signal at the finger; 

23 slope is the angle between the horizontal axis and the straight line passing through 

24 the minimum points on the left and on the right of the cardiac cycle; 

25 maxfin and minfin coincide with the systolic pressure and the diastolic pressure; 

26 and 

27 areaf is the total area of the pressure signal. 

l 13. Apparatus for measuring cardiac flow, characterized in that it comprises at 
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2 least one sensor for detecting a blood pressure signal and a computer unit 

3 connected to the said sensor for the execution of a measurement according to the 

4 method of at least one of claims 1-12, and with at least one output device of the 

5 value measured. 

1 14. Apparatus for measuring cardiac flow according to Claim 13, characterized in 

2 that the said sensor is a cardiac catheter for arterial invasive measurements. 

1 15. Apparatus for measuring cardiac flow according to Claim 13, characterized in 

2 that the said sensor is a cardiac catheter for invasive measurements in the 

3 pulmonary artery. 

1 16. Apparatus for measuring cardiac flow according to Claim 13, characterized in 

2 that the said sensor is a sensor for non-invasive measurements. 

1 17. Apparatus for measuring cardiac flow according to Claim 16, characterized in 

2 that the said sensor is a small cuff for non-invasive measurements, which is 

3 wrapped around a finger. 

1 18. Apparatus for measuring cardiac flow according to Claim 17, characterized in 

2 that the said sensor is protected during use by a second sterile cuff of the 

3 disposable type. 

1 19. Apparatus for measuring cardiac flow according to Claim 17 characterized in 

2 that said computer unit comprises a storage medium containing a computer 

3 program to execute a method according to at least one of claims 1 to12. 

1 20. Computer program loadable in a computer unit in order to execute a method 

2 according at least one of claims 1 to 12. 
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